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(Dated: October 2, 2015) 

We report the results of a search for neutrinoless double-beta decay in a 9.8 kg-yr exposure of 
130 Te using a bolometric detector array, CUORE-O. The characteristic detector energy resolution 
and background level in the region of interest are 5.1 ± 0.3 keVFWHM and 0.058 ± 0.004 (stat.) ± 
0.002 (syst.) counts/(keV-kg-yr), respectively. The median 90% C.L. lower-limit half-life sensitivity 
of the experiment is 2.9 x 10 24 yr and surpasses the sensitivity of previous searches. We find no 
evidence for neutrinoless double-beta decay of 130 Te and place a Bayesian lower bound on the decay 
half-life, T®J 2 > 2.7 x 10 24 yr at 90 % C.L. Combining CUORE-O data with the 19.75 kg-yr exposure 
of 130 Te from the Cuoricino experiment we obtain T®J 2 > 4.0 x 10 24 yr at 90% C.L. (Bayesian), the 
most stringent limit to date on this half-life. Using a range of nuclear matrix element estimates we 
interpret this as a limit on the effective Majorana neutrino mass, mpp < 270-760 meV. 


Neutrinoless double-beta (Ozx/3/3) decay is a hypothe¬ 
sized lepton-number-violating process [T] that has never 
been decisively observed. Its discovery would prove that 
lepton number is not a symmetry of nature, establish that 
neutrinos are Majorana fermions, possibly constrain the 
absolute neutrino mass scale, and support theories that 
leptons seeded the matter-antimatter asymmetry in the 
universe [2j. The clear potential for fundamental impact 
has motivated intense effort to search for this decay EHg. 

The Cryogenic Underground Observatory for Rare 
Events (CUORE) 0 0, now in the final stages 
of construction at Laboratori Nazionali del Gran 
Sasso (LNGS), promises to be one of the most sensitive 
upcoming 0 v/3/3 decay searches. The detector exploits 
the bolometric technique [2 0] in 5 x 5 x 5 cm 3 nat Te0 2 
crystals, whereby the tiny heat capacity attained by a 
crystal at ~10 mK results in a measurable increase of its 
temperature when it absorbs energy. The sought-after 
signature of Ozx/5/3 decay is a peak in the measured en¬ 
ergy spectrum at the transition energy (Qgg), which for 
130 Te is 2527.518 ± 0.013 keV HO). 

CUORE will consist of 19 towers containing 52 crys¬ 
tals each; CUORE-O is one such tower built using 
the low-background assembly techniques developed for 
CUORE [llj . The 52 crystals pT2j are held in an ultra- 
pure copper frame by polytetrafluoroethylene supports 
and arranged in 13 floors, with 4 crystals per floor. Each 
crystal is instrumented with a neutron-transmutation- 
doped Ge thermistor [[13] to record thermal pulses and 
a silicon heater to generate reference pulses [II] . The 
tower is deployed in Hall A of LNGS and exploits the 
cryogenic system, shielding configuration, and electron¬ 
ics from a predecessor experiment, Cuoricino [ME]. 

CUORE-O represents the state of the art for large- 
mass, low-background, ultra-low-temperature bolometer 
arrays. While also a competitive 0 vfi/3 decay search, it 
has validated the ultraclean assembly techniques and ra¬ 
diopurity of materials for the upcoming CUORE experi¬ 
ment. Technical details can be found in HU H2 HBHH!; 
we focus here on the first physics results from CUORE-O. 


* Present address: Physics Department, Princeton University, 
Princeton, NJ 08544, USA 
t Deceased. 


The data were collected in twenty month-long blocks 
called datasets during two campaigns which ran from 
March 2013 to August 2013 and from November 2013 to 
March 2015. For approximately three days at the begin¬ 
ning and end of each dataset we calibrated the detector 
by placing thoriated wires next to the outer vessel of the 
cryostat. Data collected between calibrations, denoted 
physics data, are used for the Ozx/3/3 decay search. 

Each thermistor voltage, except for one thermistor 
which we failed to wire bond, is continuously acquired at 
a rate of 125 Hz. Events are identified using a software 
trigger with a channel-dependent threshold of between 
30 keV and 120 keV. The trigger rate per bolometer is 
60 mHz (1 mHz) in calibration (physics) mode. Particle- 
induced pulses have rise (decay) times of ^0.05 s (^0.2 s), 
and have amplitudes of ~0.3 /rV/keV before amplifica¬ 
tion. We analyze a 5-s-long window consisting of 1 s 
before and 4 s after each trigger. The pre-trigger voltage 
establishes the bolometer temperature before the event; 
the pulse amplitude establishes the event energy. Every 
300 s, a stable current pulse is injected in each heater to 
generate tagged monoenergetic reference pulses. Noise 
waveforms are collected on all bolometers every 200 s. 

The analysis utilizes two pulse-filtering techniques, de¬ 
noted optimal filter (OF) and decorrelated optimal fil¬ 
ter (DOF), and two methods for thermal gain stabiliza¬ 
tion (TGS), denoted heater-TGS and calibration-TGS. 
The filters optimize energy resolution [21] by exploiting 
the distinct frequency characteristics of particle-induced 
vs. noise pulses. TGS corrects for small changes in the 
energy-to-amplitude response of the detection chain us¬ 
ing monoenergetic heater or calibration events. Both OF 
and heater-TGS were used for Cuoricino m- We devel¬ 
oped DOF to reduce correlated noise between adjacent 
crystals; such noise mainly affects the upper floors of the 
tower closest to cryostat noise sources [22] [23; ■ 

To recover data from the two bolometers with non¬ 
functioning heaters and from periods when temperature 
drifts in a bolometer exceeded the linear dynamic range 
of the heater-TGS, we developed calibration-TGS, which 
uses the 2615 keV 208 T1 calibration line. To successfully 
apply calibration-TGS to physics data, we monitor pa¬ 
rameters that can affect the bolometer response between 
calibrations (e.g., drifts in DC offset or amplifier gain). 
Where possible we employ both TGS methods, yielding 
up to four stabilized pulse-amplitude estimators for each 
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FIG. 1. Bottom: Energy spectra of physics (blue) and cali¬ 
bration (red) data; the latter is normalized relative to the for¬ 
mer at 2615 keV. The peaks are identified as: (1) e + e“ anni¬ 
hilation, (2) 214 Bi, (3) 40 K, (4) 208 T1, (5) 60 Co, and (6) 228 Ac. 
Top: Difference of best-fit reconstructed peak energy and ex¬ 
pected peak-energy for physics (blue points) and calibration 
(red) data. The blue line is the best-fit function to the physics 
peak residuals; the shaded band is its la uncertainty. 


event (OF and DOF, with heater- and calibration-TGS). 

To convert these to energy, we correlate prominent 
peaks in the stabilized-amplitude spectra collected in cal¬ 
ibration runs with gamma lines of known energy between 
511 keV and 2615 keV (Fig. [I]). We fit a quadratic func¬ 
tion with zero intercept to the peak-mean vs. known- 
energy points to determine a calibration function for each 
stabilized-amplitude estimator of each bolometer-dataset 
and apply these to the physics data. To avoid biasing the 
subsequent analysis we then blind the physics data in the 
region of interest (ROI) using a procedure [23j which pro¬ 
duces an artificial peak at Qpp- 

We select the best-performing energy estimator 
for each bolometer-dataset to optimize sensitivity to 
Ou/3/3 decay (quantified by the ratio of energy resolution 
of the 2615 keV calibration line to the physics data ex¬ 
posure). While the combination of OF with heater-TGS 
is the default choice, combinations involving DOF and 
calibration-TGS — which are more robust against low- 
frequency common-mode noise and long-term tempera¬ 
ture drifts, respectively — are selected if the improve¬ 
ment relative to the default is statistically significant. 
The fractions of exposure using OF with calibration- 
TGS, DOF with heater-TGS, and DOF with calibration- 
TGS are 21%, 12%, and 8%, respectively. These new 
techniques result in a 4% improvement in energy resolu¬ 
tion and a 12% increase in usable exposure. 

We select Cffi/3/3 decay candidates in the physics data 
according to the following conditions. First, we discard 


low-quality data (e.g., periods of cryostat instability or 
equipment malfunction), reducing the total exposure by 
7%. To allow a bolometer time to equilibrate after each 
event (pileup rejection) we require that the times since 
the previous event and until the next event on the same 
bolometer are greater than 3.1 s and 4.0 s, respectively. 
To reject noisy pulses which can contribute to back¬ 
ground we require each waveform to be consistent with 
a reference waveform, constructed for each bolometer- 
dataset from calibration data around the 2615 keV 208 T1 
peak. Six pulse-shape parameters characterize the wave¬ 
forms, and the acceptance criteria are tuned simultane¬ 
ously on prominent peaks in the physics data to maximize 
the signal sensitivity at each peak. These peaks range in 
energy between 146 keV and 2615 keV. The sensitivity is 
quantified by the ratio of signal accepted to square-root 
of background accepted, where the signal sample is drawn 
from events that populate each peak and the background 
is drawn from nearby off-peak events. The tuning uses 
50% of the data, randomly selected, and excludes the 
ROI. To reduce background from decays depositing en¬ 
ergy in multiple crystals (e.g., cds at crystal surfaces or 
multiple Compton scatters) we reject an event if another 
occurs in the tower within ±5 ms (anticoincidence). 

The selection efficiencies are evaluated with the frac¬ 
tion of data not used for tuning and averaged over all 
bolometer-datasets. The trigger efficiency is estimated 
from the fraction of heater pulses that produce an event 
trigger; we also exploit the heater events to measure 
the energy reconstruction efficiency (i.e., the probabil¬ 
ity for a monoenergetic pulse to reconstruct correctly). 
The combined trigger and reconstruction efficiency is 
(98.529 ±0.004)%. The combined efficiency of the pileup 
and pulse-shape selection, estimated from the fraction of 
2615 keV 208 T1 events in physics data that pass this se¬ 
lection, is (93.7 ± 0.7)%. The anticoincidence efficiency 
has two components: the probability for a Ov/30 decay to 
be fully contained in one crystal and the probability for 
it to survive accidental coincidences. The former, esti¬ 
mated from simulation [25], is (88.35 ± 0.09)%; the latter 
we find to be (99.64 ± 0.10)% using the 1461 keV q-ray 
from 40 K . The total selection efficiency is (81.3 ± 0.6)%. 

We use the high-statistics 2615 keV 208 T1 line in cali¬ 
bration data to establish the detector response to a mo¬ 
noenergetic deposit (lineshape) near the ROI. The data 
exhibit a slightly non-Gaussian lineshape characterized 
by a primary peak and a secondary peak whose mean 
is lower in energy by ^0.3% and whose amplitude is 
typically ~5% of the primary peak. Non-Gaussian low- 
energy structure was also observed in Cuoricino EMU. 
The origin of this structure in CUORE-O is under investi¬ 
gation. We studied several lineshapes, including double- 
and triple-Gaussian models; while the latter perform well 
at the 208 T1 line, we adopt the double-Gaussian lineshape 
as it is the simplest that reproduces the detector response 
over the broadest energy range. 

We parametrize the lineshape p for each bolometer- 
dataset (&, d ) as = p(pb,d,Vb,d,5b,di r lb,d)- For each 
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FIG. 2. Bottom: Calibration data near the 2615 keV 208 T1 
7 -ray line, integrated over all bolometer-datasets. The solid 
blue line is the projection of the UEML fit described in the 
main text. In addition to the double-Gaussian lineshape for 
each bolometer-dataset, the fit function includes terms to 
model a multiscatter Compton continuum, a ~ 30 keV Te 
X-ray escape peak, and a continuum background; these com¬ 
ponents, summed over all bolometer-datasets, are indicated 
by the blue dashed lines (a), (b), (c), and (d), respectively. 
Top: Normalized residuals of the data and the best-fit model. 

(6, d) pair, pb,d is the mean of the primary peak, 5b,d, 
is the ratio of the means of the secondary and primary 
peaks, Cb,d is the common Gaussian width of both peaks, 
and r]b,d is the fractional intensity of the secondary peak. 
We estimate these parameters with a simultaneous, un¬ 
binned extended maximum likelihood (UEML) fit to the 
2615 keV 208 T1 calibration line (Fig. [fj; the resulting 
best-fit parameters are denoted fib.d, <7b,d, 5b,d, and rjb,d- 
We next repeat this lineshape fit on a series of peaks 
of known energy between 511 keV and 2615 keV in the 
physics data (Fig. [lj. For a peak of known energy E, 
fj,b,d{E ) can vary around the expected calibrated energy 
via a single free parameter A n{E). To treat energy de¬ 
pendence of the resolution or possible differences in res¬ 
olution between calibration vs. physics data, we vary the 
<Jb,d relative to db,d via a global scaling parameter (E). 
We fix the 5b,d and rjb,d to the corresponding 5b,d and fjb,d- 
The energy residual parameters A fi(E) are plotted 
in Fig. |Tj A prominent outlier is the peak attributed 
to 60 Co double-gamma events which reconstructs at 
2507.6±0.7 keV, 1.9±0.7 keV higher than expected [58] ; 
a shift of 0.8 ± 0.3 keV was observed in Cuoricino 5Ti! . 
The single-escape peak of the 208 T1 2615 keV gamma at 
2104 keV also reconstructs higher by 0.84 ± 0.22 keV. 
Data taken with a 60 Co source confirm the double¬ 
gamma events reconstruct at higher energy, in agree¬ 
ment with our physics data. Simulations show their en¬ 
ergy deposit in a bolometer is less localized than the 
single-gamma lines studied; this may be responsible for 
the observed response. The double-escape peak of the 
208 T1 2615 keV line (E ~ 1593 keV) reconstructs within 
0.13 ± 0.30 keV of the expected value. Since e + e _ pairs 
and 0zz/3/3 decays share similar event topologies we as- 
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FIG. 3. Bottom: The best-fit model from the UEML fit 
(solid blue line) overlaid on the spectrum of Op/3/3 decay can¬ 
didates in CUORE-O (data points); the data are shown with 
Gaussian error bars. The peak at ~2507 keV is attributed to 
60 Co; the dotted black line shows the continuum background 
component of the best-fit model. Top: The normalized resid¬ 
uals of the best-fit model and the binned data. The vertical 
dot-dashed black line indicates the position of Qpp. 

sume the latter would reconstruct according to the cali¬ 
brated energy scale. 

We estimate the calibration offset at Qpp from a 
parabolic fit to the physics-peak residuals in Fig. [l] ex¬ 
cluding the 60 Co double-gamma and 208 T1 single-escape 
lines as outliers. We adopt the standard deviation of the 
parabolic-fit residuals as a systematic uncertainty. The 
result is A/j,(Qpp) = 0.05 ± 0.05(stat.) ± 0.09(syst.) keV. 
Similarly, fitting the resolution-scaling parameters with a 
linear function we find a a (Qpp) = 1.05±0.05. Using this 
c^a(Qpp), we estimate from calibration data the FWHM 
at Qpp of each bolometer-dataset in physics data. We 
quote the exposure-weighted harmonic mean of these 
physics FWHM values, 5.1 ± 0.3 keV, as a characteris¬ 
tic value of the detector resolution in the ROI [53]. The 
RMS of the calibration FWHM values is 2.9 keV. 

After unblinding the ROI by removing the artificial 
peak, we determine the yield of Oiz/3/3 decay events from 
a simultaneous UEML fit |26| in the energy region 2470- 
2570 keV (Fig. [3|. The fit components are: a posited 
signal peak at Qpp, a peak at ~2507 keV from 60 Co 
double-gammas, and a continuum background attributed 
to multiscatter Compton events from 208 T1 and surface 
decays [53] . We model both peaks using the established 
lineshape. For 0 v/3/3 decay, the Hb,d{Qpp) are fixed at 
the expected position (i.e., 87.00 keV — A n(Qpp) below 
flb,di where 87.00 keV is the nominal energy difference 
between Qpp and the 208 T1 line), the (Jb,d are fixed to be 
1.05 x db,d, the 5b,d and r]d,b are fixed to their best-fit cal¬ 
ibration values, and the Ojz/ 3/3 decay rate (r 0 „) is treated 
as a global free parameter. The 60 Co peak is treated 
in a similar way except that a global free parameter is 
added to the expected p,b,d to accomodate the anomalous 
double-gamma reconstruction. The 60 Co yield, although 
a free parameter, is constrained to follow the isotope’s 
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FIG. 4. Profile negative log-likelihood (NLL) curves for 
CUORE-O, Cuoricino |15H17j . and their combination. 

half-life [28] since it was cosmogenically produced above 
ground but is not replenished under ground at LNGS. 
Within the limited statistics the continuum background 
can be modeled with a zeroth-order polynomial; we con¬ 
sider first- and second-order alternatives later. 

The ROI contains 233 candidates in a total exposure of 
35.2 kg-yr of TeC> 2 , or 9.8 kg-yr of 130 Te considering the 
natural isotopic abundance of 34.167% [30] • The best- 
fit r 0 „ is 0.01 ± 0.12 (stat.) ± 0.01 (syst.) x 10 -24 yr -1 , 
and the best-fit background index in the ROI is 
0.058 ± 0.004 (stat.) ± 0.002 (syst.) counts/(keV-kg-yr). 

We evaluate the goodness of fit by comparing the value 
of the binned y 2 in Fig. [ 3 ] (43.9 for 46 degrees of free¬ 
dom) with the distribution from a large set of pseudo¬ 
experiments with 233 Poisson-distributed events in each, 
and generated with the best-fit values of all parameters; 
90% of trials return y 2 > 43.9. The data are also com¬ 
patible with this set of pseudo-experiments according to 
the Kolmogorov-Smirnov metric. We quantify the signif¬ 
icance of each of the positive and negative fluctuations 
about the best-fit function by comparing the likelihood 
of our best-fit model to the likelihood from an UEML fit 
where the fluctuation is modeled with a signal peak. For 
one degree of freedom, the most negative (positive) fluc¬ 
tuation has a probability of 0.5% (3%). The probability 
to realize the largest observed fluctuation anywhere in 
the 100-keV ROI is ~ 10%. 

We find no evidence for (ffi/3/3 decay and set a 90% 
C.L. Bayesian upper limit at r 0 „ < 0.25 x 1CT 24 yr" 1 , 
or > 2.7 x 10 24 yr (statistical uncertainties only); 
the prior used was uniform (7r(ro„) = 1 for >= 0). 
The median 90% C.L. lower-limit sensitivity for is 
2.9 x 10 24 yr. The probability to obtain a more stringent 
limit than the one reported above is 54.7%. Including 
systematic uncertainties (Table ITl) the 90% C.L. limits 
are r 0 „ < 0.25 x 1CT 24 yr" 1 or T°J 2 > 2.7 x 10 24 yr. 

To estimate systematic uncertainties we perform a 
large number of pseudo-experiments with zero and non¬ 
zero signals. We find the bias on To from the UEML 
analysis is negligible. To estimate the systematic error 
of the lineshape choice we repeat the analysis of each 
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FIG. 5. Constraints on mpp vs. lightest neutrino mass 
(miightest). For the inverted (IH, green) and normal (NH, red) 
hierarchies the central dark band is derived from the best- 
fit neutrino oscillation parameters, the lighter outer band in¬ 
cludes their 3a uncertainties [39] • The horizontal bands delin¬ 
eated by the long-dashed black lines (a), the dashed beige lines 
(b), and the dot-dashed blue lines (c) are the range of 90% C.L 
upper limits on mpp coming from (a) 130 Te (CUORE-O com¬ 
bined with Cuoricino), (b) 136 Xe (EXO-200 [40], KamLAND- 
Zen !4T] independently), and (c) 76 Ce (combined limit from 
Cerda, IGEX, HDM ;32]). The vertical arrows aim to em¬ 
phasize the range currently probed with each isotope. The 
horizontal, hashed grey band indicates the range of limits 
on mpp expected from CUORE assuming its target 90%C.L 
lower limit half-life sensitivity of 9.5 x 10 25 yr is attained. 

pseudo-experiment with single- and triple-Gaussian mod¬ 
els and study the deviation of the best-fit decay rate 
from the posited decay rate as a function of the latter. 
Similarly, we propagate the 5% uncertainty on a a (Qpp), 
the 0.09 keV energy scale uncertainty, and the choice of 
zeroth-, first-, or second-order polynomial for the back¬ 
ground. 

TABLE I. Systematic uncertainties on Foi, for zero signal (Ad¬ 
ditive) and as a percentage of nonzero signal (Scaling). 


Additive (10 24 yr 4 ) Scaling (%) 

Lineshape 

0.004 

1.3 

Energy resolution 

0.006 

2.6 

Fit bias 

0.006 

0.15 

Energy scale 

0.006 

0.4 

Bkg function 

0.004 

0.7 

Selection efficiency 


0.7% 


We combine our data with a 19.75 kg-yr exposure of 
130 Te from Cuoricino 52. The exposure-weighted mean 
and RMS FWHM energy resolution of the detectors were 
6.9 keV and 2.9 keV, respectively; the ROI background 
index was 0.169 ± 0.006 counts/(keV-kg-yr). We report 
the profile likelihoods in Fig. [4] The combined Bayesian 
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90% C.L. limit is T®J 2 > 4.0 x 10 24 yr, which is the most 
stringent limit to date on this quantity. For comparison, 
the 90% C.L. frequentist limits [3Tj are T ;| °U > 2.8 x 
10 24 yr for CUORE-O only, and Ty , 2 > 4.1 x 10 24 yr for 
the combination with Cuoricino. 

We interpret our Bayesian combined limit in the con¬ 
text of models for 0iz/3/3 decay mediated by light Ma- 
jorana neutrino exchange using the phase-space factors 
from [32], the most recent nuclear matrix element (NME) 
calculations for a broad range of models [331487] , and as¬ 
suming gA — 1-27 for the axial coupling constant. The 
resulting range for the 90% C.L. upper limit on the ef¬ 
fective Majorana mass is mg /3 < 270 650 meV; for ease 
of comparison with limits from other isotopes in the field 
(Fig-© this range excludes Ref. [35]. Including the latter 
NME, the range extends to mgg < 270 760 meV. 

In summary, CUORE-O finds no evidence for 0u/3[3 
decay of 130 Te and, when combined with Cuoricino, 
achieves the most stringent limit to date on this pro¬ 
cess. Benefiting from lower background, improved energy 
resolution, and higher data-taking efficiency, CUORE-O 
surpassed the sensitivity of Cuoricino in half the runtime. 
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